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We construct the hadronic multiplicity distribution in terms of an elementary distribution (at 
a given impact parameter) and the inelastic overlap function characterized by the observed BEL 
(Blacker-Edgier-Larger) behaviour. With suitable parametrizations for the elementary quantities, 
based on some geometrical arguments and the most recent data on e + e~ annihilation, an excellent 
description of pp and pp inelastic multiplicity distributions at the highest energies is obtained. 



I. INTRODUCTION 



^ ' It is well known that the underlying theory of hadronic interactions, the quantum chromodynamics (QCD) is not yet 
CN . able, by its own, to describe the bulk of experimental data associated with soft (long range) processess, in particular 
elastic scattering. At present, phenomenological approaches are very important as a source of information for adequate 
theoretical developments. On the other hand, from the experimental point of view, a renewed interest on both pp and 
pp elastic and inelastic scattering is expected with the advent of the next accelerator generation, the RHIC and LHC 
jjj . At this stage, due to our limited theoretical understanding of clastic scattering and, on the contrary, the success 
• of QCD in treating hard (inelastic) processes, it may be important to investigate possible connections between elastic 
and inelastic channels, even from a phenomenological point of view. In this work we shall treat this subject in the 
contexts of the impact parameter picture, unitarity and the eikonal approximation. 

Following other authors Q , we express the "complex" (overall) hadron-p multiplicity distributions (inelastic channel) 
Oh, in terms of an "elementary" distribution associated with an elementary process taking place at given impact parameter 
and the inelastic overlap function (which is constructed from elastic channel data). The novel aspects concern: (a) 
uantitative correlation between the violations of the KNO scaling j| (inelastic channel) and geometrical scaling 
(elastic channel); (b) introduction of suitable parametrizations for the elementary quantities, based on some 
geometrical arguments and the most recent data on e + e~ annihilation. With this formalism, the hadronic multiplicity 
distribution may be evaluated without any free parameter and an excellent reproduction of the experimental data on 
pp (ISR) and pp (SppS) inelastic multiplicities is achieved. A detailed discussion on all the results presented here and 
also a complete list of references to the experimental data may be found in a recent paper || . 



II. IMPACT PARAMETER FORMALISM 



In the geometrical picture, unitarity correlates the elastic scattering amplitude in the the impact parameter b space, 
r(6, s), with the inelastic overlap function, Gj„(6,s), by 2ReT(b, s) — \T(b, s)\ 2 + Gi n (b,s). For a purely imaginary 
clastic amplitude in momentum transfer space the profile function T(b, s) is real and in the eikonal approximation is 
expressed by T(b, s) = 1 — exp[— Q,(b, s)], so that 

G in {b, s) = 1 - exp[-2 tl(b, s)} = a in (b, s) (1) 

is the probability for an inelastic event to take place at b and s, namely, o~i n (s) = J d 2 h Gi n (b,s). In this picture 
the topological cross section for producting an even number N of charged particles at y/s may be expressed by 
<?n(s) = J d 2 h a N {b,s) = J d 2 h a in (b,s)a N (b,s)/a in (b,s). 

In this context, the formal connection between hadronic and elementary multiplicity distributions is obtained as 
follows. Let tp be the elementary multiplicity distribution, < n > (b, s) the average number of particles produced at b 
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and s and z — N(s)/ < n > (b, s) a KNO variable associated with the elementary process taking place at b (and s). 
Then, in general, <p =< n > (6, s)<jN(b, s)/ai n (b, s) — ip(z, s). Representing the overall multiplicity distribution by $ 
and the corresponding KNO variable by Z = N(s)/ < N > (s), where < N > (s) is the average multiplicity at y/s, 
we have in general $ =< N > (s)<7/v-(s)/(7j ra (s) = &(Z, s). Both distributions are normalized by the usual conditions 
I 

pOO nOO pOO nOO 

/ ip(z)dz = 2= / ip(z)zdz, / <5>(Z)dZ = 2 = I <$>(Z)ZdZ. (2) 
Jo Jo Jo Jo 



Now, introducting a multiplicity function as the ratio 

< n > (6, s) _ 



m(b,s), (3) 



< N > (s) 

the relationship between $ and tp follows from the above equations: 

r^^IM) ( z ) 

Jd*bG in (b,s) 9[Z ' S) - W 

This result means that, once one has parametrizations for Gi„(&, s) and the elementary quantities ip (multiplicity 
distribution) and m(b, s) (multiplicity function) the overall hadronic multiplicity distribution may be evaluated. In 
this work we consider Gi n (b,s) from analyses of elastic pp and pp scattering data (taking account of geometrical 
scaling violation) and infer the elementary quantities based on geometric arguments and experimental data on e + e~ 
annihilation, as explained in what follows. 

In the elastic channel, the breaking of Geometrical scaling is quite well described by the BEL behaviour, analytically 
expressed by the Short Range Expansion Gi n (b,s) = P(s)exp{—b 2 / AB{s)}k(x, s), with k expanded in terms of 
a short-range variable x — b exp{— ('yb) 2 / 4B (s)} . With suitable parametrizations for P(s) and B(s) an excellent 
agreement with experimental data on pp and pp elastic scattering is achieved and we shall use this well known result 
M as the input from the elastic channel. 



III. ELEMENTARY HADRONIC PROCESS 



We now turn to the discussion of the elementary hadronic processes, characterized by ip and m in Eq. (4). At a given 
impact parameter, the most elementary process is an e + e~ collision, which is a process occuring in a unique angular 
momentum state and therefore also at a given impact parameter (zero in this case). Although e + e~ annihilations 
cannot be exactly the same as collisions between hadron constituents, it is reasonable, even from a geometric point 
of view, to think that some characteristics of both processes may be similar. The point is to find out or infer what 
they could be. 

First, let us consider that the average multiplicity at given impact parameter depends on the center-of-mass energy 
in the form of a general power law < n > (b fixed, s) cx E^ M - Now, from Eq. (1), exp{— 2 f2(6, s)} is the transmission 
coefficient, i.e. the probability of having no interaction at a given impact parameter, and therefore should be 
proportional to the thickness of the target, or the energy Ecm that can be deposited at b for particle production at a 
given s. Then, we can express < n > (6, s) oc S! 7 (£>, s) and comparison with Eq. (3), allows us to infer the multiplicity 
function: 

m(M)=£0 7 (M), (5) 

with £ being calculated by the normalization condition of the overall multiplicity distribution, Eq. (2). With this, 
Eq. (4) becomes 

J d 2 h G in (b,s) 



where 



Jdb 2 G m (b, s) 
J db 2 G in (b,s)W(b,s) 



(7) 
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We proceed with the determination of the elementary distribution ip and the power coefficient 7, through quanti- 
tative analyses of e + e~ data and under the following arguments. Because the elementary process occurs at a given 
impact parameter, its elementary structure suggests that it should scale in the KNO sense. Now, since experimental 
information on e + e~ multiplicity distributions shows agreement with this scaling, we shall base our parametrization 
for ip just on these data. In particular, it is sufficient to assume a gamma distribution, normalized according to Eq. 
(2), ip(z) = 2K K z K ~ 1 exp{-Kz} /T(K). Fit to the most recent data, covering the interval 22.0 GeV < ^/s < 161 
GeV furnished K = 10.775 ± 0.064 with x7AW=508/195. 

Finally, we consider fits to the e + e~ average multiplicity through a general power law < n > e + e -— a(y / s) 7 . We 
collected the most recent experimental data at the highest energies, covering the interval 5.1 GeV < ^fs < 183 
GeV. Fitting to this equation yields a = 2.06 ± 0.02 and 7 = 0.522 ± 0.002, with x 2 /A r Dir=354/45. The above 
parametrization deviates from the data above yfs ~ 100 GeV and this contributes to the high \ 2 value. However, 
as commented before, we do not expect that e + e~ annihilation exactly represent the collisions between hadrons 
constituents. 



IV. RESULTS AND CONCLUSIONS 



With the above results we are now able to predict the hadronic multiplicity distribution &(Z, s), Eq. (6), without 
free parameters: Gi n (b, s) (and f2(6, s)) comes from analysis of the elastic scattering data and f(z) and 7 from fits to 
e + e~ data. We express $ in terms of the scaling variable Z = N / < N > where < N >= N(s) — Nq with Nq = 0.9 
leading charges removed. The predictions for pp scattering at 52.6 GeV and pp at 546 GeV are shown in Figs. 1 and 
2, together with the experimental data. The behaviour of the normalization factor £(s), as determined from Eq. (7), 
is shown in Fig. (3) in a wide interval of energy. The multiplicity functions m(b,s), Eqs. (1) and (5), corresponding 
to the energy and reactions of Figs. 1 and 2 are displayed in Fig. 4. 

From Figs. 1 and 2, we conclude that, in the context of the geometrical approach and based on the BEL behaviour 
for Gi„(b,s), our guess for the elementary hadronic quantities (based on geometrical arguments and quantitative 
analysis of e + e~ data) yields excellent predictions for the hadronic distribution <I>(Z, s). 

The fact that the average multiplicity associated with hadron constituents, increases faster then e + e~ data above 
~ 100 GeV seems quite reasonable, since at these energies we expect additional contributions from gluons/quarks 
interactions which are not present in lepton-lepton collisions. 

Anyway, despite the success of the geometrical approach in the phenomenological context, it is in general difficult to 
obtain direct connections with a microscopic theory. Some aspects of these connections have been recently discussed 
H and other results, along these lines, were also presented at this Conference M. 
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FIG. 1. Scaled multiplicity distribution for inelastic pp 
data at ISR energies compared to theoretical expectations 
using Eq. (6). 
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FIG. 2. Same as Fig. 1 for inelastic pp data at collider 
energies. 
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FIG. 3. Normalization factor from Eq. (7) as function of 
the energy. 



s" 1.0 

E 



p(bar)p - 546 GeV 

pp - 52.6 GeV 




FIG. 4. Multiplicity function from Eq. (5) for pp at 52.6 
GeV and pp at 546 GeV. 
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